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We present a study of the scaling behavior of the Renyi entanglement entropy (REE) in SU(A'^) 
spin chain Hamiltonians, in which all the spins transform under the fundamental representation. 
These SU(A'^) spin chains are known to be quantum critical and described by a well known Wess- 
Zumino-Witten (WZW) non-linear sigma model in the continuum limit. Numerical results from our 
lattice Hamiltonian are obtained using stochastic series expansion (SSE) quantum Monte Carlo for 
both closed and open boundary conditions. As expected for this ID critical system, the REE shows 
a logarithmic dependence on the subsystem size with a prefector given by the central charge of the 
SU(A) WZW model. We study in detail the sub-leading oscillatory terms in the REE under both 
periodic and open boundaries. Each oscillatory term is associated with a WZW field and decays as 
a power law with an exponent proportional to the scaling dimension of the corresponding held. We 
Hnd that the use of periodic boundaries (where oscillations are less prominent) allows for a better 
estimate of the central charge, while using open boundaries allows for a better estimate of the scaling 
dimensions. For completeness we also present numerical data on the thermal Renyi entropy, which 
equally allows for extraction of the central charge. 

I. INTRODUCTION 

Entanglement is a fascinating aspect of many-body 
quantum systems.^ It describes how a many-body wave 
function cannot in general be written as a tensor product 
of individual single-particle wave functions. The degree 
to which a wave function fails to be written as a product 
state of two subsystem wave functions can be quantified 
in terms of the entanglement entropy (EE) between these 
two subsystems. In one dimension, gapped quantum sys¬ 
tems have an EE that stays constant as the size of the 
subsystem is increased. This is consistent with the so- 
called “area law,” which states that the EE grows with 
the area of the boundary of the subsystem.^ However 
in gapless conformally invariant systems, the EE violates 
the area law and exhibits a logarithmic divergence with a 
prefactor given by the central charge.The logarithmic 
growth of entanglement has been observed in a number 
of quantum critical lattice models whose scaling limit is 
described by a GET. 

The entanglement entropy also has sub-leading terms 
that contain universal information beyond the central 
charge. Namely, one can find power law dependen¬ 
cies on the subsystem size with exponents given by 
the scaling dimensions of various conformal field theory 
(GET) operators.® Indeed, such universal contributions 
have been observed both numerically^’® and analytically.® 

These sub-leading terms in the entanglement entropy can 
be used to characterize the operator content of the un¬ 
derlying GET of various lattice models. 

In this work we consider the Renyi entanglement en¬ 
tropy (dehned in Sec. Ill) in the context of SU(A^) 

Heisenberg antiferromagnetic spin chains. We are able 
to completely characterize the underlying GET for this 
model based on the central charge from the leading log 
term, as well as the operator content contained in the 
sub-leading oscillatory terms, which decay as power laws. 

We have also measured the Renyi entropy of a subsystem 


as a function of temperature in order to gain perspective 
on finite temperature effects and as a secondary means 
of extracting the central charge. 

II. THE LATTICE HAMILTONIAN 

We consider the following Hamiltonian with the spin 
on each lattice site transforming under the fundamental 
representation of the SU(A^) algebra. 

Af-l 

= J X! (I) 

<ij> a,/3—0 

This model^® consists of a sum of bond operators that 
permute nearest neighbor spins (colors) which are labeled 
by numbers 0 through iV — 1. Hu reduces to the spin-I/2 
Heisenberg model when N = 2, and provides a natural 
extension to larger N. We will consider the case when 
J is positive (antiferromagnetic) and spins tend to anti¬ 
align with one another. Since it takes N lattice sites 
to form an SU(A^) singlet, we expect and find that the 
ground state consists of equal numbers of each color and 
is non-degenerate if the chain is an integer multiple of 
N. Finally, and most importantly for the work that we 
present here, the ground state is described by the SU(A^)i 
WZW model with central charge c = N— and N— 1 
primary fields with scaling dimensions = a — a? /N 
where a € [1, A — 

Models with SU(7V) symmetry are of both theoretical 
and experimental interest, since it has been shown that 
ultra cold atoms in optical traps can give rise to this 
symmetry.^® In fact, the model we consider here can be 
obtained from the SU(A) Hubbard model at 1/fV filling 
in the limit of large on-site repulsion (see [17,19] for a 
numerical study of the Mott transition). 

Entanglement in this class of models has previously 
been studied using DMRG for N < 4,^® though the uni- 
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FIG. 1: Comparison of QMC and exact diagonalization for 
the ground state REE of a 14 site chain with N=2 under 
periodic boundary conditions. This is shown for different 
values of the Renyi index a, as indicated in the legend. 
Colored points are data from exact diagonalization and 
black circles connected by dashed lines are obtained from 
QMC for a = 2, 3, 4, 5. Note that under periodic boundaries, 
oscillations in the REE appear when a > 1 (see the end of 
Appendix C). 


versal sub-leading oscillations were not present in the von 
Neumann entanglement entropy under closed boundary 
conditions. Other studies^^’^^ found oscillations in the 
spin-spin correlation function for different values of IV, 
verifying that the periodicity is given by 2tt jN. In fact, 
one can make a precise connection between entanglement 
and spin correlations in one dimension. 

Here we will study in detail the scaling form of these 
oscillations which are induced in the Renyi entanglement 
entropy under both open and closed boundaries. We will 
demonstrate quantitatively that the decay of these oscil¬ 
lations contain interesting information about the scaling 
dimensions of operators in the WZW field theory. 


III. NUMERICAL RESULTS 

We begin by first defining the Renyi entanglement en¬ 
tropy (REE). Take a one-dimensional system of length 
L and partition it into two segments of length I a and 
Ib = L — Ia- Construct the density matrix for the en¬ 
tire system (p), and compute the reduced density matrix 
{pa — Tr^p) by tracing over the degrees of freedom in 
Ib- The Renyi entanglement entropy is then given by 

5a(pA) = :r^log(Tr{pS}), (2) 

1 — a 

where one obtains the von Neumann entanglement en¬ 
tropy S{Ia) = ~Tr{pyi log(pA)} by taking the limit 
a —>■ 1. Appendix C reviews the extended ensemble ap¬ 
proach introduced in [24] that allows for efficient mea¬ 
surements of the REE with QMC. In Eig. 1 we com¬ 
pare our QMC results versus exact diagonalization for 
an SU(2) chain of length L = 14 under periodic bound¬ 
aries for several values of the Renyi index (a). The QMC 


delivers exact results within controllable statistical error 
bars. Similar agreement is found for N > 2 (not shown). 


A. REE periodic boundaries 


Inspired by previous work on the XXZ spin chain we 
fit our numerical data to the following scaling form® 


5„(U) = 5i°s(U) + ^°-(U) + 5„, 


where 
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where p = 1, 2 is for periodic and open boundary condi¬ 
tions, respectively. The universal parameters are the cen¬ 
tral charge c, the Fermi momentum kp, and the scaling 
dimension Ai. Fa{lA/L) is a universal scaling function 

2 Ai 

(into which a factor of | sin(fc^)|T^ has been absorbed) 
and c is a non-universal constant. In the present study 
we find that approximating Fa{lA/L) to be a constant 
allows for a sufficiently accurate extraction of the pa¬ 
rameters of interest. In the rest of this paper we take 
kp = tt/A, which can clearly be seen in the data. All 
of our simulations are performed with /3J = NL, which 
ensures that finite temperature effects are negligible. 

Fig. 2 shows our data for the REE with periodic 
boundary condition for 2 < N < 6. When N = 2, there 
is just one primary field with conformal weight Ai. In 
the case of higher A, there are primary fields with less 
relevant scaling dimensions that contribute in addition to 
the oscillatory behavior. In the periodic case, oscillations 
are very small, and one oscillatory term from the most 
relevant primary field is sufficient to describe the data. 
However, precisely because the oscillations are small, we 
are unable to reliably extract the scaling dimensions from 
our numerical fits (though it must be included for reliable 
extraction of the central charge). Higher values of the 
Renyi index indeed make it easier to extract exponents 
from the oscillations; however, this route is impractical 
since finite size effects become greater for larger a- The 
benefit of considering periodic boundaries is that it leads 
to very accurate estimates of the central charge with min¬ 
imal finite size effects, as shown in Fig. 2 and Table I. 


B. REE open boundaries, A < 4 

In order to efficiently extract the scaling dimensions 
from the REE, we consider open boundary conditions 
where oscillations are much more pronounced. We begin 
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FIG. 2: REE as a function of the subsystem size (Ia) with 
periodic boundary conditions. We set the Renyi index a = 2 
and the total chain lengths are integer multiples of N. 
Oscillations have Fermi momentum kp = n/N. Solid lines 
are best fits to the CFT scaling form. In the right inset we 
have subtracted the oscillatory and constant pieces of the 
REE, and plotted it against S]^^(Ia)/ c. This is plotted on 
top of lines with slope — 1. The best fit values for the 
central charges are given in Table I. 


N L c 

Ccft 

2 64 0.99(1) 

1 

3 74 2.01(1) 

2 

4 64 2.99(1) 

3 

5 70 3.99(1) 

4 

6 72 5.01(1) 

5 


TABLE I: Best fit central charges corresponding to Fig. 2. 
Exact values are given by ccft = — 1. These results are 

obtained by excluding the first few data points when fitting 
to the form Eq. (3). 

with N = 2,3 where there is just one distinct scaling 
dimension. The second Renyi entropy along with the 
best fit of Ai is given in Fig. 3. We find that Ai in 
the SU(2) case is not fully converged due to the presence 
of logarithmic corrections to correlations that have not 
been accounted for.^® Interestingly, this seems to have 
less of an effect in the SU(3) case where the best fit in 
the region Ia^ 1 converges close to the analytical value 
in the thermodynamic limit (see Inset of Fig. 3). 

In order to see the qualitative signature of the primary 
fields, we show in Fig. 4 the discrete Fourier transform of 
the REE appearing in Fig. 3. Before taking the Fourier 
transform we used the fact that ScHa, L) = Sc,{L—Ia, L) 
to reconstruct the REE along the entire chain length. 
We then dropped L/4 points from each edge in order to 
minimize finite size effects coming from the boundary. 
The Fourier transform is given by 

- n —1 

= ^ ( 6 ) 

V j=0 

where Sj '■ j G [0,n — 1] is the list of entries in Sa{lA) 
after the points have been dropped and n is the total 
number of elements left. We have dropped the a index 
from the discrete Fourier transform for ease of notation. 



FIG. 3: Second Renyi entropy with open boundaries for 
A = 2, 3. In the inset we plot the scaling dimension (Ai) as 
obtained by fitting the QMC data. The solid black lines are 
the exact values, and the QMC results are plotted as a 
function of the number of boundary points that are excluded 
from the fit. 
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FIG. 4: Discrete Fourier transform Eq. (6) of the REE data 
appearing in Fig. 3. Here we used the fact that 
Sa{lA,L) = Sa{L — Ia, L) to reconstruct the REE along the 
entire chain length, then dropped LjA points from each edge 
before taking the Fourier transform. Peaks in the Fourier 
transform appear at integer multiples of 2kF- 

C. REE open boundaries, A > 4 

We now move to A = 4 which is the first N for which 
there is a more than one distinct scaling dimension. We 
hence have to generalize Eq. (5) to a sum of oscillating 
terms. We use the following form, 

N-l 

SnU) = Y. f^cos{2akFlA) 

a—1 

( 7 ) 

where we achieve very good fits to our QMC data, again 
taking the universal scaling function to be a constant. 

In Figs. 5 and 6 we have fit our data from A = 4 and 
A = 5 with the oscillatory piece Eq. (7) and use it to 
extract the two distinct scaling dimensions. Eor A = 4, 5 
we find it necessary to go to even larger system sizes in 
order to show convergence of the scaling dimensions to 
their GET values. Strong finite size effects are apparent 
in the extraction of A 2 , however it is essential to include 
it as a fit parameter in order to obtain a reasonable esti- 
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FIG. 5: Second Renyi entropy with open boundaries for 
= 4 and L = 32, 64,128. The inset is similar to Fig. 3, 
although now we fit two different primary field scaling 
dimensions. Strong hnite size effects are apparent in the 
extraction of A 2 , the signature of which is much weaker 
than that of Ai. Error bars indicate stochastic error and do 
not include the systematic error inherent in the fit. 



FIG. 6: Second Renyi entropy with open boundaries for 
N = 5 and L = 30, 70,120. This figure is similar to Fig. 5. 

mate of Ai. 

Though it is difficult to obtain accurate quantitative 
estimates of A 2 using this method, clear qualitative sig¬ 
natures of all primary fields are present in the Fourier 
spectrum of the REE shown in Figs. 7 and 8. 
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FIG. 7: Fourier transform of the second Renyi entropy for 
A = 4 given in Fig. 5. 
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FIG. 8: Fourier transform of the second Renyi entropy for 
N = 5 given in Fig. 6. 

D. Effect of finite temperature on Renyi entropy 

Finally we consider the thermal Renyi entropy, which 
is defined in the same way as Eq. (2), except that the 
density matrix is no longer pure (constructed from only 
the ground state) but rather it is mixed (constructed from 
a thermal distribution of excited states): 



The thermal Renyi entropy allows us to extract the cen¬ 
tral charge at finite temperature which in practice is 
much less computationally intensive. 

Fig. 9 shows the emergence of a linear scaling region 
in the thermal Renyi entropy that now captures the en¬ 
tanglement between subsystems as well as the thermal 
entropy of subsystem A. Since the entropy goes like the 
log of the number of states, and in the high temperature 
limit all of the states are equally probable, we nat¬ 
urally expect a linear scaling region to emerge at finite 
temperature. 

We use the following scaling form to fit to our thermal 
Renyi entropy data"^’^ 

where we fix kp = tt/N. We tested that this for¬ 
mula gives the correct central charge in the linear scaling 
regime for different values of N and a and for both open 
and closed boundary conditions. Here we only present 
data for A = 3 and a = 2. 

In Fig. 10 we use the scaling form Eq. (9) to extract 
an effective central charge for A = 3 chains at different 
values of L and /3 J, which is given in the inset. We clearly 
see that the central charge flows to its analytical value in 
the limit L —>■ 00 and /3J ^ 1. In practice, we find the 
central charge approaches a value slightly higher than its 
analytical one. This is due to the fact that oscillations 
have been neglected by considering only troughs, leading 
to a linear fit with a slightly larger slope. 
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FIG. 9: Thermal Renyi entropy of a length 12 chain for 
SU(2) with both open (left) and closed (right) boundaries. 
Colored data points are obtained from exact diagonalization, 
and the black QMC data points are in perfect agreement. 
Oscillations still occur at finite temperature, but become 
drowned out at small enough values of /3J where the linear 
scaling regime emerges and open and closed chains take on 
the same scaling form. 



mation contained in the entanglement entropy. By mea¬ 
suring this quantity alone, one determines all the pa¬ 
rameters that make up the continuum description in 
terms of a CFT. One could extend this work by consider¬ 
ing sub-leading (possibly oscillating) terms in the entan¬ 
glement entropy for different representations of SU(iV). 
Such models are described by more general WZW CFTs. 
These have been studied numerically in Ref. [20] and 
Ref. [26] , though it would be interesting to see the struc¬ 
ture of oscillations that one observes and whether it is 
possible for scaling dimensions to be extracted via some 
generalization of Eq. (7). 

Partial financial support was received through NSF 
DMR-1056536. This work used the Extreme Science and 
Engineering Discovery Environment (XSEDE), which is 
supported by National Science Eoundation grant number 
ACI-1053575; in particular, resources were used on the 
Trestles cluster housed at the San Diego Supercomputing 
Center (SDSC) allocated under TG-DMR130040. Part of 
this work was completed while one author (RKK) held 
an adjunct faculty position at the TIFR. 


Appendix A: Stochastic series expansion 

In this section, we briefly review the stochastic series 
expansion (SSE) (see [27] for a comprehensive review). 
Consider a general spin Hamiltonian which is a sum of 
bond operators 

H = (Al) 

b 


FIG. 10: Thermal Renyi entropy at /3J = 15 as a function of 
subsystem size plotted for several lengths on periodic SU(3) 
chains. The inset shows effective central charges extracted 
using the form Eq. (9) for several values of total length and 
PJ. The fits only include points that are well within the 
linear scaling regime. 

IV. CONCLUSIONS 

In this paper we have investigated the Renyi entangle¬ 
ment entropy in the context of critical SU(V) spin chains 
which are described by a WZW non-linear sigma model 
in the thermodynamic limit. We showed that signatures 
of all V — 1 primary fields are present in the oscillations 
of the entanglement entropy. We further used the ana¬ 
lytical form of the oscillations given by Eq. (7) to extract 
the numerical values of the scaling dimensions, which are 
consistent with the results of CFT. 

We considered both closed and open boundary condi¬ 
tions, where the former proves effective in extracting the 
central charge, while the latter is more suitable for ex¬ 
tracting the scaling dimensions of primary fields. Finally, 
we demonstrated universal behavior of the thermal Renyi 
entropy that allows for extraction of the central charge 
with less computational effort. 

These results serve to illustrate the wealth of infor¬ 


where Hi is diagonal in the spin basis, is off-diagonal 
and b is the bond index. Here C is a constant which is 
adjusted such that the matrix elements of Hi > 0. Con¬ 
figurations in the SSE come from expanding the partition 
function, 

Z = Tr{e-^^} = 5] 5] ^{a\{Y,{Hl + Hl)r\a) 

an h 

(A2) 

on 

= (A3) 

a n {S„} 

where S'„ is a product of n diagonal and off-diagonal bond 
operators and the sum on {^n} includes all possible com¬ 
binations of such operators. We have neglected the con¬ 
stant shift (7, which should be taken into account when 
computing the energy. At this point we should mention 
that since matrix elements of diagonal operators can al¬ 
ways be made positive by the constant shift, minus signs 
only enter the SSE through off-diagonal operators. One 
is then restricted to models where H^ > 0, or models 
where the number of negative off-diagonals in any con¬ 
figuration is strictly even. The latter argument will be 
used in Appendix B to apply the SSE to our antiferro¬ 
magnetic spin Hamiltonian. 
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FIG. 11: A partition function configuration in the SSE. 
Black bars represent off-diagonal operators, and white bars 
represent diagonal operators. The operator string satisfies 
the trace condition, meaning the initial spins at the bottom 
are the same as the spins at the top. 


It is useful to introduce a pictorial representation of a 
configuration in the SSE, which is given in Fig 11. Here 
we have black (white) bars representing off-diagonal (di¬ 
agonal) operators in the spin basis, and spin colors are 
propagated through an operator string that satisfies the 
trace condition. The updating scheme for generating new 
SSE configurations will be given in the context of our 
SU(A^) model in Appendix B. 


Appendix B: QMC sampling procedure 

Previous works^^’^^ have simulated i?n by Monte- 
Carlo methods. Below we summarize how we applied 
the SSE to Hji, Eq. (1), and how MC updates of the 
configurations were carried out. 

In order to avoid minus signs associated with diagonal 
matrix elements in the SSE (see Appendix A), we shift 
the Hamiltonian by the unit matrix 

l=^|a/3)(a^| (Bl) 

Oi,j3 

h'^ = Hu- A^bondl = ^ ^(|Aa,)(ai^,| (B2) 

<ij> 

— \ail3j){aij3j\). 

We thus zero out all matrix elements between same-spin 
nearest neighbors. The off-diagonal (spin flip) opera¬ 
tors still contribute to minus signs in the SSE, thus any 
operator string that satisfies the trace and that has an 
odd number of off-diagonal operators will have a negative 
weight. On a chain with open boundary conditions, all 


configurations that satisfy the trace have an even num¬ 
ber of off-diagonal operators. This can be seen from the 
fact that there are no matrix elements between same-spin 
nearest neighbors so that each spin must return to its 
original position after propagation through the operator 
string, which requires an even number of permutations. 

With periodic boundary conditions it is possible for 
like-colored spins to trade places, which in general could 
be done by an odd number of off-diagonal matrix ele¬ 
ments. However, it can be shown that if a configuration 
has a fixed parity, i.e., the total number of each different 
type of color is either all odd, or all even, then such neg¬ 
ative weight windings are impossible. This is precisely 
what happens in the ground state of chains where the 
total length is evenly divisible by N, because the ground 
state contains L/N color singlets. Thus we will sample 
all configurations with the correct (positive) weight if we 
remain at low enough temperature to be in the ground 
state subspace and if we consider chain lengths that are 
integer multiples of N. In fact we can go further. After 
discussing the loop structure in the SSE, we will argue 
that we avoid a sign problem even at finite temperature 
in the periodic case. 

At this point we will review the Monte Carlo sampling 
procedure, highlighting the novel loop structure that al¬ 
lows for deterministic sampling of configurations. Monte 
Carlo sampling occurs in four stages, the first being a 
diagonal update, where diagonal operators are inserted 
or removed from the operator string with a probability 
given by the ratio of weights in the SSE. Next, linked 
lines are constructed that establish a loop structure to 
the configuration. A random color and starting position 
is then chosen, and one follows along the loop chang¬ 
ing colors according to the starting color until the loop 
closes. Painting loops causes the conversion between di¬ 
agonal and off-diagonal operators, and when the loop 
closes we generate a new configuration contained in the 
SSE of the partition function. Finally, one makes mea¬ 
surements on the newly generated configuration; in the 
case of the REE, checking if a transition can be made 
between different partition function ensembles (see Ap¬ 
pendix C). Loop moves are constructed so that matrix el¬ 
ements with zero weight are never generated. In our case, 
the zero matrix elements are between nearest neighbors 
with the same color. Figure 12 shows the possible vertex 
updating moves, or equivalently, how a loop update can 
change matrix elements within the SSE of the partition 
function. 

The upper boxes show updates where the incoming 
color poses no threat of a zero matrix element. The lower 
boxes require a switch and reverse loop move along with 
a loop color switch so as to avoid creating a zero matrix 
element. Once a random starting color and starting po¬ 
sition for the loop is chosen, then the path of the loop as 
well as the sequence of color changes are uniquely deter¬ 
mined by these rules. 

As we increase the value of /3J, the operator string 
grows in length (vertical direction of Fig 11), and the 
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FIG. 12: Possible vertices of the Hamiltonian and different 
loop updating moves. The left side of each box shows the 
vertex before the update, along with path and color of the 
loop (arrows). The right side of each box shows the vertex 
after it has been updated. The upper boxes show the 
“continue straight” and “switch and continue” moves where 
the incoming color will not generate a zero matrix element. 
The lower boxes are both “switch and reverse” moves that 
are required to avoid generating zero matrix elements. 

loops increase their spatial (horizontal) extent. Con¬ 
versely, loops at high temperature become localized in 
the horizontal spatial direction, and are less likely to wind 
around in the horizontal direction. It is precisely for this 
reason that the algorithm does not suffer from a sign 
problem in the case of periodic boundaries at finite tem¬ 
perature. Negative weight configurations (though possi¬ 
ble), have a negligibly small weight. As the temperature 
is lowered and the loops begin to proliferate, we begin to 
enter the total spin zero sector of the configuration space 
(the same number of every color is present) and negative 
winding configurations cannot occur. This is why it is 
crucial that we have chain lengths which are integer mul¬ 
tiples Af, since free spins (not paired in a singlet) are free 
to wander and create negative weight windings. 


Appendix C: REE and the replica trick 

We will now review the recipe for measuring the Renyi 
entanglement entropy within the context of the SSE. The 
two basic ingredients are the replica trick,^®’^® and an ex¬ 
tended ensemble QMC approach introduced in [24] . The 
Renyi entanglement entropy is given by 

5a(pA) = :r^log(Tr{pS}), (Cl) 

I — a 


where 


PA = (C2) 

is the reduced density matrix, given by tracing over the 
basis states in the B subsystem. When considering the 
entanglement between A and B at zero temperature, we 
need to have /3 sufficiently large so as to project out only 


the ground state contribution to Eq. (C2). We can ex¬ 
press the entanglement entropy in terms of a ratio of 
partition functions^^ 


^aiPA) 



log 



(C3) 


where Z'^ consists of a copies of the regular partition 
function, and has a modified trace condition in the 
A sub-region, which extends over all copies. An example 
of each type of configuration is given in Fig. 13 for the 
case a = 2. 



FIG. 13: Gonfiguration of (left) and (right). A 
transition from left to right can be made when the A spins 
in both copies match each other. Transitions from right to 
left are made when the A spins in the middle match the A 
spins in the outer trace. 

To compute the Monte Carlo average of the ratio of 
partition functions, one preforms importance sampling 
of the extended ensemble of configurations contained in 
U . Transitions between ensembles can be made 
when the spin matching condition is satisfied in the A 
sub-region. In this way, the ratio is given by the number 
of Monte Carlo steps preformed in Z^^\ divided by the 
number of Monte Carlo steps preformed in Z°‘ 



/ MCstepiZ^A^)\ 

\ Lc 


(C4) 


Sampling the extended ensemble in the SSE is partic¬ 
ularly straightforward, since transitions between ensem¬ 
bles occurs with probability 1 whenever the spin match¬ 
ing condition is satisfied in the A sub-region. This is 
due to the fact that the two configurations have identi¬ 
cal weights (number of bond operators), though this idea 
can be extended to configurations that have continuous 
degrees of freedom.In practice, the A spin matching 
condition becomes increasingly difficult to satisfy as the 
A sub-region grows, so we employ an increment trick used 
in both [29] and [24]. We note that the modified parti- 
tion function Z)^ ^ locally breaks translational invariance 
















































at a point joining the A and B subsystems.® This leads 
to oscillations in the Renyi entanglement entropy under 
periodic boundary conditions. The von Neumann entan¬ 
glement entropy lacks these oscillations under periodic 
boundary conditions as it is free from this defect. Both 


the Renyi and von Neumann entanglement entropies can 
have oscillatory terms under open boundaries due to the 
breaking of translational invariance at the edges of the 
chain. 
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